Stable spinning optical solitons in three dimensions 
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' We introduce spatiotemporal spinning solitons (vortex tori) of the three-dimensional nonlinear 

t-H ' Schrodinger equation with focusing cubic and defocusing quintic nonlinearities. The first ever found 

completely stable spatiotemporal vortex solitons are demonstrated. A general conclusion is that 
stable spinning solitons are possible as a result of competition between focusing and defocusing 

' nonlinearities. 

Oh : 

Optical solitons (spatial, temporal, or spatiotemporal) are self-trapped light beams or pulses that are supported 
by a balance between diffraction and/or dispersion and various nonlinearities. They are ubiquitous objects in optical 
media [Q. Spatiotemporal solitons (STS) ||, alias superspikes & || or light bullets were found in many works 
^J. Although they cannot be stable in the uniform Kerr (x ) medium Jfj), stability can be achieved in saturable 
quadratically nonlinear (x' 2 -*) [|[ |j, and graded-index Kerr media jo|T While a fully localized STS in three 
dimensions (3D) has not yet been found in an experiment, 2D ones were observed in a bulk x (2) medium The 
interplay of spatio-temporal coupling and nonlinearity may also play an important role in self-defocusing media [ pT[ 
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Spinning (vortex) solitons are also possible in optical media. Starting with the works 12 1, both delocalized ("dark") 
' and localized ("bright") optical vortices in 2D were investigated |L3[ [u], [l5|. In the 3D case they take the shape of a 
torus ("doughnut") fl6| , p~7|| . However, the only previously known physical model which could support stable 3D vortex 
CN| ' solitons is the Skyrme model ]l8| ] , which has recently found a new important application to Bose- Einstein condensates 
(BEC) jlj|. Our objective in this paper is to identify fundamental models of the nonlinear-Schrodinger (NLS) type 
in 3D that give rise to stable spinning solitons, as NLS models are much simpler and closer to more experimental 
• | | situations, having applications to optics, BEC, plasmas, etc. (see below). 

For bright vortex solitons stability is a major issue as, unlike their zero-spin counterparts, the spinning solitons are 
prone to destabilization by azimuthal perturbations. In 2D models with x^ and saturable nonlinearities an azimuthal 
instability was revealed by simulations |L4| and observed experimentally |l5|]. As a result, a soliton with spin 1 splits 
1 into two or three fragments, each being a moving zero-spin soliton. Simulations of the 3D spinning STS in the x^ 
model also demonstrates its instability-induced splitting into separating zero-spin solitons [|17|. Nevertheless, the x^ 
nonlinearity acting in combination with the self-defocusing Kerr (y£^) nonlinearity, gives rise to the first examples of 
stable spinning (ring-shaped) 2D solitons with spin s = 1 and 2 [j20f . It should be stressed that all the 2D spinning 
solitons actually represent static spatial beams; on the contrary, 3D solitons are moving spatiofemjwro^ ones, which 
are localized not only in the transverse plane, but also in the propagation coordinate, see below. 

A model which may support stable spinning solitons in 3D is the one with a cubic-quintic (CQ) nonlinearity, which 
(in terms of optics) assumes a nonlinear correction to the medium's refractive index in the form Sn = n-2l — n^I 2 , 
I being the light intensity. The CQ nonlinearity was derived, starting from the Maxwell-Bloch equations, for light 
propagation combining resonant interaction with two- level atoms and dipole interactions between the atoms pi), 
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or the Kerr nonlinearity of a waveguide |2^] (see also Ref. |p3[). A unifying feature of those media is competition 
between different nonlinearities, and stable solitons may exist in the range of intensities where the competition takes 
place (which may be controlled, for instance, through the density of two-level atoms in a waveguide). In fact, the 
NLS equation of the CQ type is a generic model, which also applies to Langmuir waves in plasmas [M and BEC |25f| 
(although in the latter case, three-body interactions, which give rise to the quintic term, may induce losses through 
recombination of the colliding atoms, thus making the quintic coefficient complex). 

In the first simulations of 2D solitons with spin 1 in the CQ model, it was found that they propagated in a stable 
way, provided that their energy is not too small f2"6fl . A later analysis, based on the computation of linear-stability 
eigenvalues, demonstrated that some of the spinning 2D solitons considered in Ref. Q are subject to a weak azimuthal 
instability. Nonetheless, in another part of their existence region, with very large energies, solitons with spin s = 1 
and s = 2 were confirmed to be truly stable in the 2D CQ model|2^| (all the solitons with s > 3 are unstable). 

It was recently shown by direct simulations of the CQ model p8|that 3D spinning solitons with moderate energies 
are unstable against azimuthal perturbations, while the ones with very large energies, i.e., broad "doughnuts" with a 
small hole in the center, were robust under propagation. However, a consistent stability analysis makes it necessary to 
compute eigenvalues of small perturbations. We will conclude that sufficiently broad STS with spin s — 1 are stable, 
the stability region occupying w 20% of their existence region, while all the STS with s > 2 are unstable. 

The evolution of the electromagnetic field envelope A in the dispersive CQ medium is governed by the NLS equation, 

2iK A z + V\A + k„DA tt + 2K 2 (n 2 /n )\A\ 2 A - 2K 2 {n 4 /n )\A\ 4 A = 0, t = t - z/V, (I) 

where z and t are the propagation coordinate and time, kq and V are the propagation constant and group velocity of 
the carrier wave, D > is the temporal dispersion, acting on the transverse coordinates x and y. Equation (|l|) 
does not include higher-order effects, such as self-steepening, stimulated Raman scattering, non-paraxial diffraction, 
and third-order temporal dispersion (which, in another context, were taken into regard for spatiotemporal superspikes 
in Refs. ||), as we anticipate that only broad solitons (with the temporal width ~ I ps), for which these effects are 
small, may be stable. 

Defining rcscalcd variables u — y/n^Jn^A, T = n 2 \J 2ko /Dnon4T, Z = (ko^I/^o^) z , and (X, Y) = 
K o n 2-\/2/^o n 4 ( x i y)i we transform Eq. ([j]) into a normalized form p8[ , 

iuz + (u X x + uyy + u T t) + \u\ 2 u - \u\ A u = 0. (2) 

STS solutions to Eq. (Q) are sought as u — U(r, T) exp(is8) exp(inZ), where r and 9 are the polar coordinates in the 
transverse plane, k is a propagation constant parameterizing the family of solutions sought for, and s is an integer 
spin. The real amplitude U obeys the equation 

(U rr + r^Ur - s 2 r- 2 U + U TT ) -kU + U 3 - U 5 =0, (3) 

supplemented by the condition that U must decay exponentially as r — > oo and T — > oo (due to the definition of t 
in Eq. (0), the latter condition implies that the solution's snapshot taken at t = const is localized in the propagation 
coordinate z - in fact, exactly the same way as ID solitons are localized in optical fibers |29||). 
Equation (||) conserves the energy E = J J j \u(X,Y,T)\ 2 dXdYdT, Hamiltonian 

H = JJ J [\ux\ 2 + \u Y \ 2 + \uT\ 2 -(l/2)\u\ 4 + (l/3)\u\ e }dXdYdT, (4) 

momentum (equal to zero for the solutions considered), and angular momentum in the transverse plane, L = 
III (d(f)/d8) \u\ 2 dXdYdT , where 4> is the phase of the complex field u. Relations between L, H and E for a 
stationary spinning soliton follow from Eq. (§): L = sE; H = kE - % f J 2irrU 6 (r, T)drdT @. 

We have numerically found families of 3D spinning solitons with a toroidal shape. To quantify the solutions, in Fig. 
I we show the propagation constant k and H vs. E for both s = and s = 1, 2 solitons. They exist for E exceeding 
a threshold value, which increases with s. The full and dashed lines in Fig. 1 correspond to stable and unstable 
branches according to results presented below. The s = branch of the solutions is divided into stable and unstable 
portions on the basis of the known criterion which states that the fundamental (s = 0) soliton branch undergoes a 
change in the stability where dE/dn = p0[ . However, this criterion ignores azimuthal instability, which is frequently 
fatal for spinning solitons. 

The most revealing information on the stability of solitons is provided by analysis of a linearized version of Eq. (|^) . 
To this end, we seek perturbation eigenmodes of the general form, 

u{Z,r,T,9) -U(r,T)exp[i{s6 + kZ)} = f(r, T) exp {X n Z + i[(s + n)6 + kZ]} 

+g*(r, T) exp {X* n Z + i[(s - n)0 + kZ}} , (5) 
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where n > is an arbitrary integer azimuthal index of the perturbation, A„ is the (complex) instability growth rate 
sought for, and the functions / and g obey equations 

B 2 f B 2 f Bf 

*W + Qf~ 2 + q£ + r" 1 ^ - (s + nfr- 2 f - nf + (2 - W 2 ) U 2 f + (l - 2U 2 ) U 2 g = 0, 
-iK9 + + ^) 2 ^ 2 5 - Kg + (2 - 3t/ 2 ) f/ 2 5 + (1 - 2C/ 2 ) U 2 f = 0. (6) 



r 



The solutions must decay exponentially at r — > oo, and vanish as rl s± ™ at r — > 0. 

To solve Eqs. (||), we used a known numerical procedure Q, [H], which produces results presented in Fig. 2. 
The most persistent unstable eigenmode has n = 2, for both s = 1 and s = 2. As is seen in Fig. 2, with the 
increase of k, the instability of the soliton with s = 1, accounted for by ReA2, disappears at k = K st « 0.13, and 
the stability region extends up to k = K^^l t ~ 0.17, corresponding to infinitely broad solitons (which implies that 
the vortex of the dark-soliton type ]12] |, that may be regarded as an infinitely broad spinning soliton, is stable too). 

The relative width of the stability region is ( K^] t — K st ^ /^ofSt ~ However, there is no stability region for 3D 

solitons with s — 2, in contrast to the 2D vortex solitons in the CQ model [^Tj. In the case when a spinning soliton 
is unstable, its instability is oscillatory; the corresponding frequency, ImA, is of the same order of magnitude as ReA 
at the maximum-instability point (see Fig. 2), and A becomes purely imaginary at k = K s t- 

The above results were checked in direct simulations of Eq. (|J) by means of the Crank-Nicholson scheme combined 
with the Gauss-Seidel iteration procedure. In Fig. 3 we show the amplitude and energy vs. Z for the soliton with 
s=l, generated by two different initial configurations, with the same energy of the initial configuration, Eq — 13070. 
Robustness of the spinning STS is attested to by the fact that it can be generated from a Gaussian with a nested 
vortex whose shape is far from the soliton's exact form. (Energy loss evident in Fig. 3(b) is caused by emission 
of radiation in the course of the formation of the stable STS; naturally, the loss is larger for the initial Gaussian 
configuration, which is farther from the exact soliton's shape.) Figure 4 shows the gray-scale contour plots of the 
intensity and phase of both the input Gaussian with a nested vortex and emerging spinning STS at Z = 400. 

The instability of the s = 2 solitons is illustrated in Figs. 5 and 6. The azimuthal instability breaks them into 
zero-spin solitons which fly out tangentially, relative to the circular crest of the original soliton. It is noteworthy 
that, at an early stage of the evolution shown in Fig. 5, the spinning soliton splits into two fragments, in accordance 
with the growth-rate calculations predicting the dominant instability to be against the n = 2 perturbation mode, 
but the subsequent nonlinear evolution is more involved. In Fig. 6(b), the spinning STS splits into three fragments 
with unequal energies. The fragmentation of this soliton (having k — 0.09) into three parts is in accordance with the 
growth-rate calculation, see Fig. 2(b). 

In conclusion, we have found the first example of stable three-dimensional spinning solitons in a dispersive medium 
which combines cubic and quintic nonlinearities. Only sufficiently broad solitons with spin s = 1 may be stable. How- 
ever, the existence of stable spinning 3D solitons is a generic fact, as it is not limited to the cubic-quintic nonlincarity: 
our preliminary studies indicate that these stable physical objects may also occur in media with competing quadratic 
and self-defocusing cubic nonlinearities. In fact, a condition which turns out to be necessary for the existence of stable 
spinning solitons is a competition between two different nonlinearities, one focusing and the other one defocusing. 
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FIG. 1: The propagation constant k (a) and Hamiltonian H (b) of the 3D spinning soliton vs. its energy E. 

FIG. 2: The growth rate of perturbations, ReA, with different values of the azimuthal number n (indicated by labels near the 
curves) vs. the soliton's propagation constant k: (a) s = 1; (b) s = 2. 
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FIG. 3: Evolution of the maximum amplitude of the soliton with s = 1 (a) and its energy (b), as generated by initial 
configurations in the form of a Gaussian with a nested vortex (continuous curves) or a torus close to the stationary spinning 
soliton corresponding toK = 0.15 (dashed curves). 
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FIG. 4: The formation of the soliton with spin s = 1: (a) the initial Gaussian with a nested vortex; (b) its phase field; (c) the 
spinning soliton at Z = 400; (d) the phase field at Z — 400. A cross section of the fields at T — is shown. 



FIG. 5: Gray-scale plots showing the developing instability of the spinning soliton with s — 2 and k = 0.13 at Z = (a), 
Z = 600 (b), Z = 620 (c), and Z = 640 (d). 



FIG. 6: Isosurface plots illustrating the instability of the s — 2 soliton with k = 0.09: (a) Z = 0, (b) Z = 250. 



